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Public transport systems in 22 Polish cities have been analyzed. Sizes of these networks range 
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I. INTRODUCTION 

Since the explosion of the complex network science 
that has taken place after works of Watts and Stro- 
gatz pj as well as Barabasi and Albert 0, II] a lot of 
real- world networks have been examined. The exam- 
ples arc technological networks (Internet, phone calls 
network), biological systems (food webs, metabolic sys- 
tems) or social networks (co-authorship, citation net- 
works) 0, IE IE G}- Despite this, at the beginning 
little attention has been paid to transportation net- 
works - mediums as much important and also sharing 
as much complex structure as those previously listed. 
However, during the last few years several public trans- 
port systems (PTS) have been investigated using vari- 
ous concepts of statistical physics of complex networks 
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Chronogically the first works regarding transportation 
networks have dealt with power grids QJ 0, IE 13 O nc 
can argue that transformators and transmission lines 
have little in common with PTS (i.e. underground, 
buses and tramways), but they definitely share at least 
one common feature: embedding in a two-dimensional 
space. Research done on the electrical grid in United 
States - for Southern California 0, 

ill 

and for the 

whole country [lOj as well as on the GRTN Italian power 
network revealed a single-scale degree distributions 
(p{k) ex exp(— ah) with a rs 0.5), a small average connec- 
tivity values and relatively large average path lengths. 

All railway and underground systems appear to share 
well known small- world properties [lj- Moreover this 
kind of networks possesses several other characteristic 
features. In fact Latora and Marichiori have stud- 
ied in details a network formed by the Boston subway 
[l2l. IT3I Il4l ] . They have calculated a network efficiency 
defined as a mean value of inverse distances between 
network nodes. Although the global efficiency is quite 
large i? 9 ; & = 0.63 the local efficiency calculated in the 
subgraphs of neighbors is low Ei oca i = 0.03 what indi- 
cates a large vulnerability of this network against acci- 
dental damages. However, the last parameter increases 
to E [ocal = 0.46 if the subway network is extended by 



the existing bus routes network. Taking into account ge- 
ographical distances between different metro stations one 
can consider the network as a weighted graph and one is 
able to introduce a measure of a network cost. The es- 
timated relative cost of the Boston subway is around 0.2 
% of the total cost of fully connected network. 

Sen et al. ^E] have introduced a new topology de- 
scribing the system as a set of train lines, not stops, and 
they have discovered a clear exponential degree distribu- 
tion in Indian railway network. This system has shown 
a small negative value of assortativity coefficient. Seaton 
and Hackett |16| have compared real data from under- 
ground systems of Boston (first presented in [l4| ) and Vi- 
enna with the prediction of bipartite graph theory (here: 
graph of lines and graph of stops) using generation func- 
tion formalism. They have found a good correspondence 
regarding value of average degree, however other proper- 
ties like clustering coefficient or network size have shown 
differences of 30 to 50 percent. 

In works of Amaral, Barrat, Guimera et al. [E 
Ipgl Il9| a survey on the World-Wide Airport Network 
has been presented. The authors have proposed trun- 
cated power-law cumulative degree distribution P(k) oc 
k~ a f(k/k x ) with the exponent a — 1.0 and a model of 
preferential attachment where a new node (flight) is in- 
troduced with a probability given by a power-law or an 
exponential function of physical distance between con- 
nected nodes. However, only an introduction of geo- 
political constrains 01 (i- e - only large cities are allowed 
to establish international connections) explained the be- 
havior of betweenness as a function of node degree. 

Other works on airport networks in India .20] and 
China |2l| have stressed small-world properties of those 
systems, characterized by small average path lengths 
({/) » 2) and large clustering coefficients (c > 0.6) with 
comparison to random graph values. Degree distribu- 
tions have followed either a power-law (India) or a trun- 
cated power-law (China). In both cases an evidence of 
strong disassortative degree-degree correlation has been 
discovered and it also appears that Airport Network of 
India has a hierarchical structure expressed by a power- 
law decay of clustering coefficient with an exponent equal 
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to 1. 

In the present paper we have studied a part of data 
for PTS in 22 Polish cities and we have analyzed their 
nodes degrees, path lengths, clustering coefficients, as- 
sortativity and betweenncss. Despite large differences 
in sizes of considered networks (number of nodes ranges 
from N = 152 to N = 2881) they share several universal 
features such as degree and path length distributions, log- 
arithmic dependence of distances on nodes degrees or a 
power law decay of clustering coefficients for large nodes 
degrees. As far as we know, our results are the first com- 
parative survey of several public transport systems in the 
same country using universal tools of complex networks. 




II. THE IDEA OF SPACE L AND P 



FIG. 2: Map of examined cities in Poland. 



To analyze various properties of PTS one should start 
with a definition of a proper network topology. The idea 
of the space L and P, proposed in a general form in |l5| 
and used also in |lj| is presented at Fig. The first 
topology (space L) consists of nodes representing bus, 
tramway or underground stops and a link between two 
nodes exists if they are consecutive stops on the route. 
The node degree k in this topology is just the number 
of directions (it is usually twice the number of all PTS 
routes) one can take from a given node while the distance 
I equals to the total number of stops on the path from 
one node to another. 



plan of the city roads has been mapped into an " informa- 
tion city network" . In the last topology a road represents 
a node and an intersection between roads - an edge, so 
the network shows information handling that has to be 
performed to get oriented in the city. 

We need to stress that the spaces L and P do not take 
into account Euclidean distance between nodes. Such 
an approach is similar to the one used for description 
of several other typ es of network systems: Internet Hi , 
power grids |l~oL 1 1 1| . railway or airport networks [2(J, 



(a) 



(b) 




line A 
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FIG. 1: (Color online) Explanation of the space L (a) and the 
space P (b). 

Although nodes in the space P are the same as in the 
previous topology, here an edge between two nodes means 
that there is a direct bus, tramay or underground route 
that links them. In other words, if a route A consists 
of nodes a,, i.e. A = {a±, a,2, a n }, then in the space 
P the nearest neighbors of the node a\ are a,2, a^, a„. 
Consequently the node degree k in this topology is the 
total number of nodes reachable using a single route and 
the distance can be interpreted as a number of transfers 
(plus one) one has to take to get from one stop to another. 

Another idea of mapping a structure embedded in two- 
dimensional space into another, dimensionless topology 
has recently been used by Rosvall et al. in where a 



III. EXPLORED SYSTEMS 

We have analyzed PTS (bus and tramways systems) 
in 22 Polish cities, located in various state districts as 
it is depicted at Fig. [2J Table [I] gathers fundamental 
parameters of considered cities and data on average path 
lengths, average degrees, clustering coefficients as well as 
assortativity coefficients for corresponding networks. 

Numbers of nodes in different networks (i.e. in differ- 
ent cities) range from N = 152 to N = 2811 and they 
are roughly proportional to populations I and surfaces S 
of corresponding cities (see Fig. |3J. One should notice 
that other surveys exploring the properties of transporta- 
tion networks have usually dealt with smaller numbers of 
vertices, such as N = 76 for U-Bahn network in Vienna 
0, N = 79 for Airport Network of India (ANI) [2l[, 
N = 124 in Boston Underground Transportation System 
(MBTAlQ or N = 128 in Airport Network of China 
(ANC) pq. Only in the case of the Indian Railway Net- 
work (IRN) [HI where N = 579 and World-Wide Airport 
Network (WAN) with 3880 nodes sizes of networks 
have been similar or larger than for PTS in Poland. Very 
recently, von Ferber et al. p3| have presented a paper on 
three large PTS: Dusseldorf with TV = 1615, Berlin with 
N = 2952 and Paris where TV = 4003. 
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TABLE I: Data gathered on 22 cities in Poland. S stands for the surface occupied by the city (in km 2 ) (241 ], I is the city's 
population in thousands of inhabitants |24| and N is the number of nodes (stops) in the network. (I) is the average path 
length, (k) - the average degree value, c is the clustering coefficient and r - the assortativity coefficient. Indexes L and P stand, 
consequently, for the space L and the space P. Properties of parameters defined in spaces L and P will be discussed in sec. 
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IV. DEGREE DISTRIBUTIONS 

A. Degree distribution in the space L 

Fig. 0]shows typical plots for degree distribution in the 
space L. One can see that there is a slightly better fit to 
the linear behavior in the log-log description as compared 
to semi-logarithmic plots. Points k = 1 are very peculiar 
since they correspond to routes' ends. Remaining parts 
of degree distributions can be approximately described 
by a power law 

p(k) ~ k-' 1 (1) 

FIG. 3: Log-log plot of the dependence of number of nodes although the scaling cannot be seen very clearly and 

N on surface S (circles) and population I (triangles). it is limited to less than one decade. Pearson correlation 

coefficients of the fit to Eq. range from 0.95 to 0.99. 
Observed characteristic exponents 7 are between 2.4 and 
4.1 (see Table |TT|l, with the majority (15 out of 22) 7 > 3. 
Values of exponents 7 are significantly different from the 
value 7 = 3 which is characteristic for Barabasi- Albert 
model of evolving networks with preferential attachment 
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FIG. 4: Degree distributions in the space L for four chosen 
cities. Plots (a) and (b) show the distributions in log-log scale 
while plots (c) and (d) - in the semi-log scale. 



and one can suppose that a corresponding model for 
transport network evolution should include several other 
effects. In fact various models taking into account effects 
of fitness, atractiveness, accelerated growth and aging of 
vertices [26J or deactivation of nodes |27l |28| lead to 7 
from a wide range of values 7 £ (2, 00). One should also 
notice that networks with a characteristic exponent 7 > 4 
are considered topologically close to random graphs [2{| 
- the degree distribution is very narrow - and a difference 
between power-law and exponential behavior is very sub- 
tle (see the Southern California power grid distribution 
in presented as a power-law with 7 w 4 and in 
depicted as a single-scale cumulative distribution). 

Degree distributions obtained for airport networks are 
also power-law (ANC, ANI) or power-law with an expo- 
nential cutoff (in the case of WAN). For all those systems 
exponent 7 is in the range of 2.0 — 2.2, which differs sig- 
nificantly from considered PTS in Poland, however one 
has to notice, that airport networks are much less de- 
pendent on the two-dimensional space as it is in the case 
of PTS. This effect is also seen when analyzing average 
connectivity ((k) = 5.77 for ANI, (k) = 9.7 for WAN and 
(k) = 12 — 14 for ANC depending on the day of the week 
the data have been collected). 

Let us notice that the number of nodes of degree k = 1 
is smaller as compared to the number of nodes of degree 
k = 2 since k = 1 nodes are ends of transport routes. 
The maximal probability observed for nodes with degree 
k = 2 means that a typical stop is directly connected 
to two other stops. Still some nodes (hubs) can have a 
relatively high degree value (in some cases above 10) but 
the number of such vertices is very small. 



B. Degree distribution in the space P 

In our opinion, the key structure for the analysis of 
PTS are routes and not single bus/tramway stops. There- 
fore we especially take under consideration the degree 
distribution in the space P. 

To smooth large fluctuations, we use here the cumula- 
tive distribution P{k) Q according to the formula 



P(k) = J p(k)dk (2) 
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FIG. 5: (Color online) P(k) distribution in the space P for 
eight chosen cities 

The cumulative distributions in the space P for eight 
chosen cities are shown at Fig|SJ Using the semi-log scale 
we observe an exponential character of such distributions: 

P(k) = Ae- ak (3) 

As it is well known Q the exponential distribution © 
can occur for evolving networks when nodes are attached 
completely randomly. This suggests that a correspond- 
ing evolution of public transport in the space P pos- 
sesses an accidental character that can appear because 
of large number of factors responsible for urban develop- 
ment. However in the next sections we show that other 
network's parameters such as clustering coefficients or 
degree-degree correlations calculated for PTS are much 
larger as compared to corresponding values of randomly 
evolving networks analyzed in 3] . 

In the case of IRN UM de gree distribution in the space 
P has also maintained the single-scale character P(k) ~ 
e~ ak with the characteristic exponent a — 0.0085. Val- 
ues of average connectivity in the studies of MBTA 
((k) = 27.60) and U-Bahn in Vienna ((fc) = 20.66) are 
smaller than for considered systems in Poland, however 
one should notice that sizes of networks in MBTA and 
Vienna are also smaller. 
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TABLE II: Coefficients 7 and a with their fitting errors A7 
and Aa. Fitting to the scaling relation Q has been performed 
at whole ranges of degrees k while fitting to 10 has been 
performed at approximately half of available ranges to exclude 
large fluctuations occurring for higher degrees (See Fig. |5j. 



C. Average degree and average square degree 

Taking into account the normalization condition 
P{k m in) = 1 we get the following equations for the aver- 
age degree and the average square degree: 



(k) = 
(k 2 ) = 
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Dropping all terms proportional to e ak " 
simplified equations for (k) i (k 2 ): 
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"4 (5) 
we receive 
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Since values of k m i n range between 3 and 16 and 
they are independent from network sizes N as well as 
observed exponents a we have approximated k m i n in 
Eqs. (JSJ - (J7J) by an average value (mean arithmetical 
value) for considered networks, (k m i n ) « 8.5. At Figs. 
El and H we present a comparison between the real data 
and values calculated directly form Eqs. © and 0. 
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FIG. 6: (k) as a function of a. Circles are real data values, 
while the line corresponds to Eq. (JSJ 
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FIG. 7: (k 2 ) as a function of a. Circles are real data values, 
while the line corresponds to Eq. iJTJ 



V. PATH LENGTH'S PROPERTIES 



A. Path length's distributions 



Plots presenting path length distributions p(l) in 
2k 2 spaces L and P are shown at Figs. [S] and El respectively. 

(k 2 ) - k 2 mm + = + ^ (7) 

The data well fit to asymmetric, unimodal functions, fn 
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fact for all systems a fitting by Lavenberg - Marquardt 
method has been made using the following trial function: 



p(l) = Ale- Bf+Cl 
where A, B and C are fitting coefficients. 



(8) 
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FIG. 8: Fitted path length distribution in the space L 

Inserts at Figs. |Hl and present a comparison between 
experimental results of {I) and corresponding mean val- 
ues obtained from Eq. JHJ . One can observe a very good 
agreement between averages from Eq. (JSJ and experi- 
mental data. 

The agreement is not surprising in the case of Fig. 
since the number of fitted data points to curve © is quite 
small, but it is more prominent for Fig. [S] 
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FIG. 9: Fitted path length distribution in the space P. 

Ranges of distances in the space L are much broader as 
compared to corresponding ranges in the space P what is 
a natural effect of topology differences. It follows that the 
average distance in the space P is much smaller ((/) < 3) 
than in the space L. The characteristic length 3 in the 



space P means that in order to travel between two differ- 
ent points one needs in average no more than two trans- 
fers. Other PTS also share this property depending on 
the system size the following results have been obtained: 
(/) = 1.81 (MBTA), (l) = 1.86 (Vienna), (l) = 2.16 
(IRN). In the case of the space L the network MBTA 
with its average shortest path length (I) = 15.55 is plac- 
ing itself among the values acquired for PTS in Poland. 
Average path length in airport networks is very small: 
(I) = 2.07 for ANC, (I) = 2.26 for ANI and (/) = 4.37 
for WAN. However, because flights are usually direct (i.e. 
there are no stops between two cities) one sees imme- 
diately that the idea of the space L does not apply to 
airport networks - they already have an intrinsic topol- 
ogy similar to the space P. Average shortest path lengths 
(I) in those systems should be relevant to values obtained 
for other networks after a transformation to the space P. 

The shape of path length distribution can be explained 
in the following way: because transport networks tend to 
have an inhomogeneous structure, it is obvious that dis- 
tances between nodes lying on the suburban routes are 
quite large and such a behavior gives the effect of ob- 
served long tails in the distribution. On the other hand 
shortest distances between stops not belonging to subur- 
ban routes are more random and they follow the Gaussian 
distribution. A combined distribution has an asymmetric 
shape with a long tail for large paths. 

We need to stress that inter-node distances calculated 
in the space L are much smaller as compared to the num- 
ber of network nodes (see Table I). Simultaneously clus- 
tering coefficients cl are in the range (0.03,0.15). Such 
a behavior is typical for small-world networks .1] and the 
effect has been also observed in other transport networks 
U 13 HI 13 E3 The small world property is 
even more visible in the space P where average distances 
are between (1.80, 2.90} and the clustering coefficient cp 
ranges from 0.682 to 0.847 which is similar to MBTA 
(c = 0.93), Vienna (c = 0.95) or IRN (c = 0.69). 



B. Path length as function of product kikj 

In |2^| an analytical estimation of average path length 
(/) in random graphs has been found. It has been shown 
that (I) can be expressed as a function of the degree dis- 
tribution. In fact the mean value for shortest path length 
between i and j can be written as [^j : 



lij {ki ? kj ) — 



■In kikj +ln((fc 2 ) - (A;)) + lnA- 7 



ln((fc2)/(fc)-l) 



1 

2 

(9) 

where 7 = 0.5772 is Euler constant. 

Since PTS are not random graphs and large degree- 
degree correlation in such networks exist we have as- 
sumed that Eq. @ is only partially valid and we have 
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FIG. 10: (Color online) Dependence ol hj on kikj in the space 
L 



written it a more general form 



(hj) = A - B log kikj. 



(10) 



To check the validity of Eq. lfTU|l we have calculated 
values of average path length between kj as a function 
of their degree product kikj for all systems in the space 
L . The results are shown at Fig. ^3 which confirms the 
conjunction (|10fl . A similar agreement has been received 
for the majority of investigated PTS. Eq. I|1U|) can be 
justified using a simple model of random graphs and a 
generating function formalism |34| or a branching tree 
approach [U El IH, IH- In fact the scaling relation 
(|10[) can be also observed for several other real world 
networks [13,1111111111. 

It is useless to examine the relation IjlOU in the space P 
because corresponding sets kj consist usually of 3 points 
only. 



VI. CLUSTERING COEFFICIENT 

We have studied clustering coefficients Ci defined as a 
probability that two randomly chosen neighbors of node 
i possess a common link. 

The clustering coefficient of the whole network seems 
to depend weakly on parameters of the space L and of 
the space P. In the first case its behavior with regard 
to network size can be treated as fluctuations, when in 
the second one it is possible to observe a small decrease 
of c along with the networks size (see Table HJ. We shall 
discuss only properties of the clustering coefficients in the 
space P since the data in the space L are meaningless. 

It has been shown in that clustering coefficient in 
IRN in the space P decays linearly with the logarithm of 
degree for large k and is almost constant (and close to 



unity) for small k. In the considered PTS we have found 
that this dependency can be described by a power law 
(see Fig. [nj: 



c(k) 



-0 



(11) 



Such a behavior has been observed in many real sys- 
tems with hierarchical structures [H H^. In fact, one 
can expect that PTS should consist of densely connected 
modules linked by longer paths. 
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FIG. 11: c(k) for Gdansk (triangles) GOP (squares) and 
Warszawa (circles). Dashed lines are fits to Eq. 1 1 U with 
following exponents: Gdansk - (3 = 0.93 ± 0.05, GOP - 
(3 = 0.81 ± 0.02 and Warszawa - (3 = 0.57 ± 0.01. All data are 
logarithmically binned with the power of 1.25. 

Observed values of exponents (3 are in the range (5 £ 
(0.54,0.93). This can be explained using a simple ex- 
ample of a star network: suppose that the city transport 
network is a star consisting of n routes with L stops each. 
Node i, at which all n routes cross is a vertex that has the 
highest degree in the network. We do not allow any other 
crossings among those n routes in the whole system. It 
follows that the degree of node i is ki = n(L — 1) and the 
total number of links among the nearest neighbors of i 
is Ei = n(L — 1)(L — 2)/2. In other words the value of 
the clustering coefficient for the node with the maximum 
degree is: 



c(fcn 



2Ei 



L-2 



k t (h - 1) n(L - 1) - 1 



(12) 



where k max = n{L — 1). It is obvious that the minimal 
degree in the network is k m i n = L — 1 and this correspon- 
dences to the value c(k m i n ) = 1. Using these two points 
and assuming that we have a power-law behavior we can 
express (3 as: 



(3 



lnc(k max ) - lnc(fc rl 



In 



In k r . 



In k rt 



L-2 
n(L-l)- 

lnn 



(13) 



Because n(L — 1) 3> 1 and L — 1 « X — 2 we have (3 w 1. 

In real systems the value of clustering coefficient of the 
highest degree node is larger than in Eq. (|i2|l due to mul- 
tiple crossings of routes in the whole network what leads 
to a decrease of the exponent (3 (see Fig. 1111) . This de- 
crease is also connected to the presence of degree-degree 
correlations (see the next Section). 



VII. DEGREE-DEGREE CORRELATIONS 

To analyze degree-degree correlations in PTS we have 
used the assortativity coefficient r, proposed by Newman 
[35l | that corresponds to the Pearson correlation coeffi- 
cient j^Bl of the nodes degrees at the end-points of link: 



(a) 



(b) 



Si Ji^i M Si Ji Si ^» 



(14) 



where M - number of pairs of nodes (twice the number 
of edges), ji, ki - degrees of vertices at both ends of i- 
th pair and index i goes over all pairs of nodes in the 
network. 
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FIG. 12: The assortativity coefficient r in the space P as a 
function of TV 

Values of the assortativity coefficient r in the space L 
are independent of the network size and are always pos- 
itive (see Table [J, what can be explained in the follow- 
ing way: there is a little number of nodes characterized 
by high values of degrees k and they are usually linked 
among themselves. The majority of remaining links con- 
nect nodes of degree k = 2 or k = 1, because k = 2 is an 
overwhelming degree in networks. 

Similar calculations performed for the space P lead to 
completely different results (Fig. I12fl . For small networks 
the correlation parameter r is negative and it grows with 
N, becoming positive for N > 500. The dependence can 
be explained as follows: small towns are described by 
star structures and there are only a few doubled routes, 




line A 
line B 
line C 
line D 






FIG. 13: (Color online) Crossing of 4 routes of 5 stops each, 
(a) In the star example there is only one hub and assortativity 
coefficient is equal to r = —0.25 according to Eq. I|15[l . In 
case (b) a few hubs exist due to a multiple crossing of routes 
and r = —0.19. Upper diagrams - the space L, lower diagrams 
- the space P. 



so in this space a lot of links between vertices of small 
and large k exist. Using the previous example of a star 
network and taking into account that the degree of the 
central node is equal to k c = n(L — 1), the degree of any 
other node is k = L — 1, after some algebra we receive 
the following expression for the assortativity coefficient 
of such a star network: 



L - 1 



(15) 



Let us notice that the coefficient r is independent from 
the number of crossing routes and is always a negative 
number. 

On the contrary, in the large cities there are lots of con- 
nections between nodes characterized by large k (trans- 
port hubs) as well as there is a large number of routes 
crossing in more than one point (see Fig. I13|) . It follows 
that the coefficient r can be positive for such networks. 
A strange behavior for the largest network (GOP) can 
be explained as an effect of its peculiar structure: the 
system is rather a conglomerate of many towns than a 
single city. Thus, the value of r is lowered by single links 
between the subsets of this network. 

At Fig. E| we show coefficients (3 as a function of r in 
the space P. One can see that in general positive values 
of the assortativity coefficient correspond to lower values 
of (3, being an effect of existence of several links between 
hubs in the networks. 

Reported values of assortativity coefficients in other 
transport networks have been negative (r = —0.402 for 
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ANI |2g and r = -0.033 for IRN and since these 
systems are of the size N < 600 thus it is in agreement 
with our results. 




FIG. 14: j3 coefficient (see Eq. lllU 'l as a function of the 
assortativity coefficient r in the space P. 



VIII. BETWEENNESS 

The last property of PTS examined in this work is be- 
tweenness t 3.9| which is the quantity describing the "im- 
portance" of a specific node according to equation |40j ]: 



cr jk 



(16) 



where, (Xjk is a number of the shortest paths between 
nodes j and k, while (Jjk{i) is a number of these paths 
that go through the node i. 



A. Betweenness in the space L 

Fig. E| shows dependence of the average betweenness 
(g) on node degree calculated using the algorithm pro- 
posed in 0] ( see also E3)- Data at Fig. El fit well to 
the scaling relation: 
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FIG. 15: The average betweenness (g) as a function of k in 
the space L for three chosen cities. 
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FIG. 16: r\ coefficient as a function of network size N 



distribution thus it can suggest that large PTS are more 
random than small ones. Such an interpretation can be 
also received from the Table |n] where larger values of the 
exponent 7 are observed for large cities. 



B. Betweenness in the space P 



9 ~ k» (17) 

observed in Internet Autonomous Systems 0], co- 
authorship networks 0] and BA model or Erdos-Renyi 
random graphs |4Cj . 

The coefficient r\ is plotted at Fig. El as a function 
of network size. One can see, that 77 is getting closer 
to 2 for large networks. Since it has been shown that 
there is r\ = 2 for random graphs 40] with Poisson degree 



The betweenness as a function of node degree k in the 
space P is shown at Fig. El O ne can see large differ- 
ences between Fig. 1151 andHTI In the space P there is a 
saturation of (g) for small k what is a result of existence 
of the suburban routes while the scale-free behavior oc- 
curs only for larger k. The saturation value observed in 
the limit of small k is given by {g(k m i n )) = 2(N — 1) and 
the length of the saturation line increases with the mean 
value of a single route's length observed in a city. 
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FIG. 17: The average betweenness (g) as a function of k in 
the space P for three chosen cities. 



IX. CONCLUSIONS 



of transfers, the degree distribution is exponential with 
characteristic exponents a = 0.013 — 0.050. Distributions 
of paths in both topologies are approximately given by a 
function p(l) = Ale~ Bl +cl . Small world behavior is ob- 
served in both topologies but it is much more pronounced 
in space P where the hierarchical structure of network is 
also deduced from the behavior of c(k). The assortativ- 
ity coefficient measured in the space L remains positive 
for the whole range of N while in the space P it changes 
from negative values for small networks to positive values 
for large systems. In the space L distances between two 
stops are linear functions of the logarithm of their degree 
products. 

Many of our results are similar to features observed 
in other works regarding transportation networks: un- 
derground, rail way or airline systems 0, 0, 0, 0, 0, 
13 El El El Hi EJ El- All such networks tend to 
share small-world properties and show strong degree- 
degree correlations that reveal complex nature of those 
structures. 



In this study we have collected and analyzed data for 
public transport networks in 22 cities that make over 
25 % of population in Poland. Sizes of these networks 
range from N = 152 to N ~ 2881. Using the concept 
of different network topologies we show that in the space 
L, where distances are measured in numbers of passed 
bus/tramway stops, the degree distributions are approx- 
imately given by a power laws with 7 = 2.4 — 4.1 while 
in the space P, where distances are measured in numbers 
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